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Abstract 

A special class of composite laminates composed of soft rubbery matrices and stiff 
reinforcements made of steel wires or synthetic fibers is examined, where each constituent 
behaves in a nonlinear fashion even in the small strain domain. Composite laminates made 
of plies stacked at alternating small orientation angles with respect to the applied axial 
strain are primarily dominated by the nonlinear behavior of the reinforcing fibers. How- 
ever; composites with large ply orientations or those perpendicular to the loading axis, 
will approximate the behavior of the matrix phase and respond in even a more complex 
fashion for arbitrarily stacked plies. The geometric nonlinearity due to small cord rotations 
during loading was deemed here to have a second order effect and consequently dropped 
from any consideration. The user subroutine USRMAT within the Micromechanics Analy- 
sis Code with the Generalized Method of Cells (MAC/GMC), was utilized to introduce 
the constituent material nonlinear behavior. Stress-strain behavior at the macro level was 
experimentally generated for single and multi ply composites comprised of continuous 
Nylon-66 reinforcements embedded in a carbon black loaded rubbery matrix. Compari- 
sons between the predicted macro composite behavior and experimental results are excel- 
lent when material nonlinearity is included in the analysis. In this paper, a brief review of 
GMC is provided, along with a description of the nonlinear behavior of the constituents 
and associated constituent constitutive relations, and the improved macro (or composite) 
behavior predictions are documented and illustrated. 


1.0 Introduction 

In flexible or soft composites, the high modulus, low elongation cords carry most of 
the load, and the low modulus, high elongation rubber matrix preserves the integrity of the 
composites and transfers the load. The primary objective of this type of composite is to 
sustain large deformation and fatigue loading while providing high load carrying capacity. 
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The classical approach of establishing the constitutive law for such composites is based on 
the traditional and well developed anisotropic theory of stiff laminates for small linear 
elastic deformations. Many reviews have demonstrated the suitability of this approach in 
predicting the elastic properties and the overall behavior of these composites under differ- 
ent types of deformation fields 1 . In this paper, the effort is focused on capturing the mate- 
rial nonlinearity by an alternative micromechanical modeling of composites approach; that 
is, one which utilizes the mechanical and physical characteristics of both constituents to 
describe the macro behavior of the composite. 

Recently, the Method of Cells (MOC) 2,3,4 and its extension the Generalized Method 

of Cells (GMC) 5,6,7 have emerged as an efficient micromechanical approach that can pre- 
dict the micro and macro performance of composites. The method accounts in an explicit 
fashion for the geometric and mechanical attributes of the constituents, a feature that dis- 
tinguishes GMC from all other averaging schemes (micromechanical models). These 
capabilities are captured in software packages such as the Micromechanics Analysis Code 

(MAC/GMC) 8 , and the First -Order Nonlinear Elasticity (FONE) 2 code. These codes are 
structured in a modular fashion that allow a user to implement additional constitutive 
models through a user defined subroutine. It is this additional feature in MAC/GMC that 
the present work capitalizes upon through the implementation of a nonlinear elastic con- 
stitutive model for each constituent. 

NASA Glenn Research Center has developed a user-friendly software package (MAC/ 
GMC) based on the generalized method of cells (GMC) and coupled it with classical lam- 
ination theory. The basic formulation of this micromechanical model expresses in explicit 
fashion the influence of fiber (cord) diameter, inter fiber distance, intra ply gauge, and the 
constituents’ mechanical properties on the global (macro) and local (micro) behavior of 
the composites. The code was tested against diverse composite configurations and results 

were compared to available experimental results with favorable agreement 7,9 . More 
recently, it was found that the linear elastic model included in MAC/GMC was adequate 
in performing analysis of elastomeric composites reinforced with steel wires with the 
applied load in the wire direction. However, for synthetic cords such as nylon, polyester, 
and rayon, and for cases where the load is orthogonal or oriented at a large angle with 
respect to the direction of the reinforcing cords, the predicted macro results using linear 
elastic assumptions deviated noticeably from the experimental data. This is due primarily 
to the nonlinear behavior of the reinforcing cords and the rubbery matrix. To circumvent 
this discrepancy, a nonlinear elastic material model for the reinforcing cord and two sepa- 
rate constitutive laws for the rubber matrix were formulated by accounting for incom- 
pressibility and the observed softening behavior; all nonlinear models where implemented 
into MAC/GMC via its user defined subroutine USRMAT. 

The present paper begins with a brief mathematical description of GMC in section 2. 
In section 3 the relevant constituent constitutive laws are described. Finally, in section 4, 
using both linear and nonlinear constitutive models, the numerical predictions obtained 
from MAC/GMC and their corroboration with experimental data are documented. 
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2.0 Generalized Method of Cells (GMC) 


GMC can predict the elastic and inelastic thermomechanical response of both continu- 
ous and discontinuous multiphased composite materials with an arbitrary internal micro- 
structure and reinforcement shape. It is a continuum-based micromechanics model that 
provides closed-form expressions for the macroscopic composite response in terms of the 
properties, size, shape, distribution, and response of the individual constituents or phases 
that make up the material. These constituent materials can be represented using any elastic 
and/or viscoplastic deformation and life (e.g., continuum damage mechanics fatigue 
model) model. The periodic nature of composites typically allows one to identify a repeat- 
ing unit cell that can be used as a building block to construct the entire composite. The 
properties of this unit cell are thus representative of the properties of the entire assem- 
blage. In the original MOC the unit cell consisted of a single fiber (cord) subcell sur- 
rounded by three matrix subcells, however in GMC this has been generalized such that the 
unit cell is subdivided into an arbitrary number of subcell phases. In Fig. 1 a periodic 
assemblage and repeating unit cell are illustrated, from which the approximate description 
of the composite is developed. 





where (3, y identify the particular subcell 



3 


Figure 1. RYE for a continuous fiber reinforced composite. 
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The basic approach has been presented elsewhere by Aboudi and co-workers (4 * 8) ; 
here only the basic steps required are briefly outlined and have been simplified to address 
linear elastic material behavior and mechanical loading only. The relations derived are for 
unidirectional fiber reinforced materials with pseudo square fibers extending in the x ( 
direction and arranged in a doubly periodic array in the x 2 and x 3 directions (see Fig. 1). 
Given the repeating unit cell, identified in Fig. 1 by dotted lines, and the four new local 
_((3) _(y) 

coordinate systems ( , .*3 ) introduced at the center of each ( 3 , y subcell, one begins 

by constructing a first order theory in which the displacements are expanded linearly in 
terms of the distances from the center each subcell. That is, 

,(Py) _ JPY) ^ ,(P)JPY) , *(Y)y(PY) J = l 23 (]) 


u i 


= w. 


+ x 2 9 . + * 3 


where y^fPY) are the displacement components at the center of each subcell, and, 
(PY) (Py) 

9 / and \j / • characterize the linear dependence of the displacements on the local 

coordinates, (see Fig. 1 ). 


The components of the small strain tensor are given as follows: 

p (PY) _ l r -) „(PY), a f .(PY)i • • , 

ij ~2 [6 <j d Ji ] 1,J= ,2,3 


( 2 ) 


where 3, = 


dx* 


d, = 


dx 
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Next the classical linear elastic constitutive relationship is assumed for both fiber and 
matrix so that the corresponding stiffness matrix is strain independent, that is. 


.(Py) _ c (Py) £ (Py) 


( 3 ) 


(py) 

where C is the elastic stiffness matrix expressed in terms of the engineering con- 
stants of the subcell. 


The average stresses in the composite, O • • , are determined from the average stresses, 




, in the subcells, i.e.; 
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where 

and 


Vr 


v = (h.+h.y 


, which is the area of the repeating unit cell, 
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To establish the composite constitutive relations, the continuity of tractions and dis- 
placements are imposed at the boundaries, in an average sense, between the constituents. 
The continuity of the displacements at the interfaces of each subcell will provide relation- 
ships between the average stresses and average strains in the composites and the microva- 
riables in the subcells of the representative cell. The conditions for the continuity of 
tractions provide the average normal stresses after lengthy mathematical manipulations. 
For the sake of completeness and yet without the loss of simplicity, only a brief examina- 
tion of the governing equations is provided. 


The condition of continuous normal and tangential displacements at the interfaces of 
the RYE are imposed on an average basis between subcells (1 y) and (2y), for example: 


I 


u 


(iy) 


dx 


(Y) 


- 0 ) 
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u 


(2y) 


(2) _ h 2 
Xl ~ 2 


« 3 (Y) 


( 6 ) 


From these conditions, the average strain in the subcells can be obtained. The average 
strains in the composite are given by: ~ 




\_ 

V 



v e (pY) 

Py £ 0' 


(7) 


P,Y= 1 


This equation establishes relationships between the average stresses and microvariables in 
the subcells of the RVE. 

The condition of equal average normal, axial and transverse stresses across the inter- 
faces provides needed algebraic equations to determine all unknown microvariable param- 
eters. Eventually, the following composite constitutive relations are established in the 
form: 


G = Be 


( 8 ) 
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where 0 are the average macro composite stresses, £ are the uniform applied macro com- 
posite strains, and B is the effective elastic stiffness tensor defined in terms of the concen- 


(PY) (0y) 

tration matrix A of the micro subcell, and the elastic stiffness tensor C of each 

subcell: 


(PyMPy) 


b = -L- y y h.h y c^ v A 

h l h 2 ^ ^ P Y 


(9) 


P = 1 Y = 1 


This was a brief summary of the equations that form the basis of the micro-to-macrome- 
chanics analysis which describe the behavior of a heterogeneous media. A more generalized 
version of this method consisting of, discontinuous as well as continuous fiber reinforce- 
ment subcells, inelastic material response, and thermomechanical load histories, has been 

fully developed and implemented within MAC/GMC 8 . MAC/GMC has uniquely enhanced 
the basic capabilities of GMC in numerous ways. The one which will be exploited in this 
work, is the fact that a variety of constituent constitutive models can be accessed via a 
library of constituent constitutive models or conveniently implemented through a user defin- 
able subroutine. The following is a demonstration of how one can modify the linear elastic 
model provided in MAC/GMC into a more elaborate nonlinear elastic model that can better 
describe the mechanical response of the constituents. 


3.0 Characterization of Nonlinear Elastic Constitutive Laws 

Currently MAC/GMC provides, in addition to three viscoplastic models (Bodner-Par- 
tom, Robinson Viscoplastic, and Generalized Viscoplasticity with Potential Structure) and 
one linear transversely isotropic model, a user-definable constitutive model that can be 
applied to either (or both) fiber and matrix constituents in a composite. Subsequently, when- 
ever the use of a linear (strain-independent) model is described, this implies use of the built 
in transversely isotropic model supplied within MAC/GMC, yet specialized to its isotropic 
equivalency. To deduce the fundamental constants such as the axial and shear moduli, a least 
square approximation of the experimental data for each constituent was employed; wherein 
the modulus of elasticity of the rubber compound was determined to be 1.87 ksi and that of 
the reinforcing cord to be 282.36 ksi. Consequently the strength of anisotropy is approxi- 
mately 151. The specification and characterization of the nonlinear (strain-dependent) 
reversible models for both cord and rubber matrix are described below. 

3.1 Fiber Reinforcement 

The composite samples used for this study, were manufactured with multi filament 
Nylon 66 cords twisted in 3 strands with a 5x5 twist per inch (each strand is twisted with 5 
twists per inch in a direction that gives it an S shape appearance, then the cord is formed by 
twisting all 3 strands with 5 twists per inch in an opposite direction that gives it a Z shape 
appearance) and a matrix made of a moderately carbon black loaded natural rubber. The 
axial behavior of the reinforcing cord was determined from the experimental macroscopic 
stress-strain response under strain control with a strain rate of 12 in/in/min; see Fig. 2. 
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Strain 


Figure 2. Nonlinear behavior of a 1260/3 Nylon cord 

In the absence of an explicit constitutive law for the nylon cord, a phenomenological 
nonlinear model was utilized; wherein the cord was idealized as a structure in and of itself. 
The nonlinearity was introduced by considering the resulting discretized macroscale cord 
modulus to be strain dependent as shown in Table 1. The corresponding nonlinear fit is 
also shown in Fig. 2. 


TABLE 1. Cord Moduli as function of the strain level 


Strain level £ 

Modulus-Msi 

0.0 - 0.008 

- 17.25e + 0.700 

0.008-0.012 

- 15.50e + 0.686 

0.012-0.016 

- 15.758 + 0.689 

0.016-0.02 

- 14.258 + 0.665 

0.02 - 0.024 

- 1 1.758 + 0.615 

0.024 - 0.0.7 

- 5.508 + 0.465 
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TABLE 1. Cord Moduli as function of the strain level 


Strain level £ 

Modulus-Msi 

0.03 - 0.04 

-3.70e + 0.411 

0.04 - 0.05 

- 1.30e + 0.315 

0.05 - 0.06 

0.10e + 0.245 

0.06 - 0.08 

1.206 + 0.179 

0.08-0.10 

1.756 + 0.135 


Note, this is only a partial description of how the cord modulus changes over the range of 
the applied strain. One can extend the description over the expected domain of application. 


3.2 Rubber Matrix 

Rubber compounds are widely used as the continuous phase in elastomeric compos- 
ites. Although, their behavior is constrained by the presence of the reinforcing fibers 
(cords), there are instances where the overall composite behavior is predominantly dic- 
tated by that of the rubber, e.g., when the composite is loaded transverse to the fiber direc- 
tion. Essentially, these materials can undergo large deformations while remaining elastic. 
The continuum mechanics theory used to model such behavior is called hyperelasticity. 
This theory can be formulated in terms of principal invariants and principal stretches with 
an additional parameter to treat the incompressibility of the material. The Mooney-Rivlin 

and Ogden models* l0) are best suited for very large deformations. Alternatively, in the 
small deformation regime, the classical isotropic linear elasticity model, expressed with a 
bulk modulus and a shear modulus as a function of strain, is a very convenient approach to 
implement in MAC/GMC. 

3.2.1 Classical Elasticity Extended 

Isothermal linear elasticity assumes a linear relationship between the Cauchy stress 
and the total infinitesimal strain. But, with proper modification, the nonlinear behavior of 
the rubber is captured by establishing the dependency of its elastic stiffness tensor on the 
level of strain magnitudes. In addition, the incompressibility of the rubber was captured in 
the hydrostatic component of the total stress by assigning a value close to 0.5 for the Pois- 
son’s ratio. The following traditional stress decomposition into two components is uti- 
lized: 


°u= «o> 

where p is the hydrostatic stress component and is related to the volumetric strain £^. 
and the bulk modulus, K, by the following relations: 
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p = Kz with K = 


3(1 -V) 


on 


where E represents the Young’s modulus and V the Poisson’s ratio of the given material. 

The deviatoric stress, S (; -, components are related to the associated deviatoric strain 
components e-j and the shear modulus G as follows: 

S iJ = C€y U2) 

wherein, the deviatoric strains are related to the total strains, Z - , and the volumetric 
strains according to the following relationship: 


e U ~ E ij~ £ kk b ij 


( 13 ) 


For computational purposes, a more useful form of the constitutive law can be pre- 
sented in tensorial form as follows: 


® ij C ijkl £ kl 


( 14 ) 


where the elastic stiffness tensor is presented in matrix form as follows: 


K + ^G K-\g K~\g 0 0 0 
K-\g K + ^G K-^G 0 0 0 

K -^G K-^G K + ^G 0 0 0 

0 0 0 G 0 0 

0 0 0 0 G 0 

0 0 0 0 0 G 


(15) 

The nonlinearity of the rubbery material was added to this elastic Hookean constitutive 
law by expressing the shear modulus as a function of the applied strain magnitude, so as to 
account for the network decay and compound breakdown. The strain dependency of the 

modulus of the carbon black filled rubber was first reported by Payne 1 1 but it was also the- 
oretically approached, by Kraus 12 . Herein, the elastic modulus G is described using an 
Arrhenius type relationship characterizing its exponential decay. 
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(16) 


G = 


G 0 e 


(-ay) 


where G 0 is the initial modulus in the absence of any previous history, y is the applied 

strain amplitude, and Ot is a phenomenological constant that describes the dependency of 
the modulus of the material on the strain softening and the irreversible breakdown under 
static strain. G 0 and Ot are material parameters that can be determined experimentally. 


To generalize the above uniaxial dependency of the shear modulus to that of a multiax- 
ial deformation field, eq. (16) can be rewritten as: 

~ ., (-<*,/(/,- 3)) 

G = G Q e (17) 

where /, = (X, + A, + A 3 ) ■ Now given an incompressible or nearly incompressible body 

2 | 

(G = A|A.->A-, = 1 ) under pure shear deformation (i.e., A, = A , A? = 1, A^ = — ) it can 

'A 2 

be shown that / j becomes, 


, ,2 , 1 
/ 1 - ^ + 1 + — 


(18) 


Also, given the fact that the strain in the torsional deformation mode considered in the 
analysis can be approximated in terms of the stretch ratio A as. 


y ~ tany 



(19) 


such that upon squaring the stretch, i.e., y = / j - 3 , one can clearly see that eq. (17) 
reduces to the uniaxial expression given in eq. (16). 

In order to characterize eq. (16), a laboratory testing procedure was implemented 
using systematic sweeps of increasing strain amplitude performed on parallel plate speci- 
men subjected to torsional shear. A characterization temperature of 25°C and a frequency 
of 1 Hz were selected on the Rheometrics System IV to achieve the best reproducibility 
over the entire strain range. The correlation delay was set to 10 seconds, so that a 3 to 5 
cycle break-in would occur prior to the response measurement. From the 10 experiments, 
275 stress strain data pairs were obtained, from which a subset of data representing shear 
modulus verses strain was established (see Figure 3). Such a set of data is usually suffi- 
cient to perform a non-linear least squares regression, leading to the determination of the 

parameters G 0 and a . A numerical method based on the Marquardt-Levenberg approach 
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was used to perform non-linear regression to obtain the two parameters G 0 and (X . For a 
typical cured rubber compound at room temperature these two constants are Gq= 5.237 

ksi, and OC = 14.7 1 3. 



- ■ ■ - 

0.00 0.10 0.20 0.30 0.40 0.50 

Strain 


Figure 3. Shear Moduli as a function of strain level 
3.2.2 Mooney-Rivlin 

To demonstrate the versatility of the software package MAC/GMC, and to provide an 
alternative way to implement a given constitutive law, a strain energy based model was 
introduced and compared against the previous model. The polynomial Mooney-Rivlin 
strain energy function was chosen to account for the material nonlinearity and incompress- 
ibility of the rubbery matrix as follows: 

W = C, (/j - 3) + C 2 (/ 2 - 3) + K(I^ - 3) (20) 

where C b C 2 are two material constants, K represents the material incompressibility, and 
/ j J 2'^ ^ are t ^ iree stra i n invariants defined as follows: 

/, = trie) = c — 

1 2 = j(tr 2 (c)-tr(c 2 )) = \ic kk c kk -c jk c kj ) 
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h = de '\ c ij\ <20 

where c is the Lagrangian strain tensor defined in terms of the deformation tensors 
T 

( c = F F ) and expressed as function of the strain components as follows: 


£ 




( 22 ) 


The associated Second Piola-Kirchhoff stress tensor (S-) is obtained by differentiat- 
ing with respect to the strain using the chain rule as follows: 



dw_ = (dW^J_ dw^_ 2 dw ^_ 3 > 

dc ij l a/ l 9c /; + dl 2 dc ij + dI 3 dc ijj 


(23) 


by substituting for the derivatives of the strain invariants and simplifying the following 
expression is obtained: 

S ij = 2 ( c l 8 y + C 2^\ b ij ~ c ip + def\c.j\ (C (24) 

Two hypothetical cases of an RVE containing four subcells, with properties taken to be 
those of the matrix, were analyzed with the two types of material models described above. 
The same strain-softening function for the matrix material was implemented into the 
USRMAT formulation. The macro stress-strain calculated by MAC/GMC is plotted in 
Fig. 4 and illustrates the excellent agreement between both nonlinear models and the 
experimental data. 
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Figure 4. Comparison between the classical and strain energy material models 


4.0 Results 

By incorporating in MAC/GMC either type of nonlinear constitutive law for the rub- 
ber behavior (including softening and incompressibility features) in addition to the strain 
dependency of the cord modulus the most accurate composite response is expected. To 
demonstrate this, an experimental program was conducted for single and multi-ply lami- 
nates. Comparison of the experimental results with numerical predictions of the composite 
behavior, based on the previously discussed constituent fiber/matrix constitutive laws, is 
discussed below. Note, all macro calculations represent true predictions as only constitu- 
ent cord and matrix behavior were correlated with experimental data and a strong bond 
was assumed to be present between cord and matrix phases. 

4.1 Unidirectional Composite 

A set of single ply unidirectional specimens with the dimensions L 12” x W 1” x T 
0.035” were built with each having approximately a 24 percent fiber volume ratio and a 
preferred direction of either [0], [30], [45] or [90] degrees. A gray square was drawn on 
the surface of each sample in the middle region where strains were measured in the princi- 
pal directions. Displacements and induced axial and lateral forces were recorded by a 
video camera and a load cell respectively. The load-displacement data from the composite 
samples were converted into stress-strain data for a direct comparison with the output gen- 
erated by MAC/GMC. A MAC/GMC input file 8 is listed in Appendix A for the case 
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where the mechanical properties of the reinforcing cord were assumed linear, i.e., to be 
constant and independent of strain magnitude. 

The 12 inch sample length was adopted after conducting an aspect ratio study* in 
which little to no change in the uniformity of the strain field at the center of the sample 
was observed for specimens of this length and longer. In addition, the various correction 

factors where examined for the elastic properties discussed in the literature* l4) for off-axis 
laminae testing. It was determined that the specimen aspect ratio of length and width uti- 
lized was sufficient so that no correction associated with axial-shear coupling was 
required. 

Furthermore, a study was conducted to investigate the influence of the assumption of 
periodicity which is at the core of the MOC and GMC formulation. Consequently, com- 
posite laminates with the same cord volume fraction (V=0.247) were cured into one, two 
and four [0] ply laminates. The samples were cut into similar dimensions as discussed 
above and tested under the same loading conditions. The macro stress-strain responses 
were recorded and compared with little to no difference being observed; due to the large 
modulus mismatch between cord and matrix and the fact that the response of the lamina is 
dominated by the cord. 

The samples were uniaxially stretched with an Instron Model 1122. The lateral and 
longitudinal strains were measured by an image processing system which include a CCD 
video camera, a NuVision image processing unit and software developed by Perceptics 
Corporation. An optical target was drawn at the center of the samples with a silver pen as 
shown in Fig. 5. The target was deformed into a lozenge shape when the specimen was 
stretched. The maximum and minimum dimensions of the lozenge, which can be con- 
verted into longitudinal and lateral strains, were obtained during the measurement. Each 
test was repeated 3 times with a new specimen, the repeatability was reasonably good, and 
the average values were reported. 



Figure 5 Test fixture 
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Comparison between the experimental and numerically predicted composite tensile 
response for all four cord orientations, i.e., [0], [90], [30] and [45] are shown in Figs. 6-8. 
Two predictions are shown per orientation, one employing a linear constituent constitutive 
model and the other a nonlinear (strain-dependent) for both cord and matrix phase. As one 
might expect, the importance of including nonlinearity in a given phase (cord and/or 
matrix) is highly dependent upon the ply orientation. Considering, the [0] degree laminate, 
as shown in Fig. 6, it is clear that in the small strain domain between 0.0 and 0.02, the 
maximum discrepancy of the linear model is on the order of 4.0%, and increases to a max- 
imum of 1 2.0% for the larger strain magnitude of 0. 1 . This discrepancy, however, is signif- 
icantly reduced by accounting for the nonlinearity of the cord behavior as in this case the 
composite response is dominated by that of the cord. Alternatively, for the case when the 
applied strain direction is perpendicular (i.e., transverse) to the cord direction the rubbery 
matrix phase now becomes the dominate constituent influencing the composite response, 
see Fig. 7. Examining Fig. 7, it is clear that utilizing either a classical or Mooney-Rivlin 
nonlinear model (implemented through the user defined material model) for the rubbery 
matrix provides the most accurate composite tensile prediction. Also, the assumption of a 
strong fiber/matrix bond, as stipulated in the analysis, is clearly justified since the compos- 
ite response is stiffer (although only slightly due to the low volume fraction) than that of 
the matrix alone. 

The additional feature related to the rubber incompressibility was not fully explored to 
the same depth as the softening effect. It is believed that under the extensional deformation 
mode investigated in this paper, the primary source of nonlinearity is affecting the shear 
modulus of elasticity more than the bulk modulus. Therefore, the results reported in Fig. 7 
were generated with a Poisson’s ratio of 0.49 and a constant bulk modulus of K= 10.475 
ksi. 



Figure 6. Stress-strain comparison for [0] degree ply. 
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Similar improvements in the accuracy (when using a nonlinear material model for the 
constitutents) of the predicted composite behavior as compared to that measured, are 
observed for the two off-axis (i.e., [30] and [45]) tests with respect to applied strain direc- 
tion, shown in Fig. 8. In the case of the [30] laminate, the linear cord and matrix model 
will start deviating from the actual data at approximately 0.04 strain, whereas when the 
nonlinear cord and matrix model are utilized the range of accuracy is increased up to 
approximately 0.08 strain with the discrepancy between linear and nonlinear predictions 
verses observed data being a factor of two at 0. 12 strain, where the nonlinear model under- 
predicts the experimental data by approximately 12% and the linear over-predicts by 
approximately 26%. A similar trend (although not shown in Fig. 8) was noted for the other 
angles with the discrepancy between linear and nonlinear model predictions versus exper- 
imental data for the [90] ply orientation being as much as a factor of six. 

Such lack of accuracy at large strain values is expected, since for large orientation 
angles the amount of cord reorientation becomes significant and the overall response is 
dominated not only by the nonlinear behavior of the cord and matrix but also the geomet- 
ric nonlinearity and its effects on the stiffness of the composites. 



Figure 7. Stress-strain comparison for [90] degree ply. 
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Figure 8. Stress-strain comparison for off-axis orientations. 

Since lateral as well as longitudinal strains were recorded, a comparison between the 
measured and predicted effective (composite) Poisson’s ratio can be made and is shown in 
Figs. 9 and 10. Predictions were made using both the generalized nonlinear classical elas- 
ticity and the Mooney-Rivlin constitutive model for the matrix, described in section 3.2. 
Clearly, the predictive ability of MAC/GMC is very good for all angles of orientation, 
with a maximum deviation from experimental results of approximately 5 percent. Further- 
more, although both cord and matrix constituent properties possess differing degrees of 
strain dependent nonlinearities; only in the case of [0] and [90] laminae is the effective 
Poisson’s ratio a function of applied strain when using the classical elasticity model. Also, 
Figs. 9 and 10 demonstrate the fact that only when one examines multiple stress or strain 
components (here the e 2 2 and e ll are represented by the effective Poisson’s ratio) can a 
distinction be made between multiaxial representations. For example, in Figs. 9 and 10 we 
see that the classical elasticity generalization appears to represent the experimental data 
better than does the Mooney-Rivlin approach discussed in section 3.2.2, as employing the 
Mooney-Rivlin constitutive model produces a strain dependent effective Poisson’s ratio 
for all ply orientation angles. This is in contrast to the observed experimental behavior. 


N ASA/TM— 1 999-209066 


17 





Applied Strain 


Figure 9 Comparison of predicted and measured effective Poisson’s ratio for [0], 
[30] and [45] ply orientations. 
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Figure 10 Comparison of predicted and measured effective Poisson’s ratio for 
[90] ply orientation. 
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4.2 Cross Ply Laminate 


A similar analysis was conducted for a cross ply laminate. Here the samples consisted 
of four plies with the dimensions L 5.5” x W 0.7874” x T 0.22”, stacked at 10 degree 
angles with respect to the applied strain direction in an alternating fashion, i.e, +[10], - 
[10], +[10], -[10] degree. The samples were first conditioned by subjecting them to 10 
loading-unloading cycles to overcome the history imposed on the sample by the mold dur- 
ing the cure cycle, then tested under axial strains. Two example listings of the revised 
MAC/GMC version 2.0 input file are presented in Appendix B. The first utilizes the lin- 
ear elastic model supplied within MAC/GMC to represent both the reinforcing cord and 
matrix. The material properties are supplied by the user as indicated by setting the param- 
eter IDB=n. The second illustrates how one would specify the needed parameters to 
invoke the user defined subroutine. 

The experimentally obtained load-extension data was once again converted into stress- 
strain data so as to conform with the output generated by MAC/GMC. As in the unidirec- 
tional laminate case, both linear and nonlinear constituent models were applied. As 
expected, after subjecting the cross ply laminate to the same loading condition as in the 
unidirectional case, a noticeable improvement in the predictive capability of MAC/GMC 
was realized when nonlinear constituent descriptions were employed, as shown in Fig. 1 1 . 
The nonlinear cord description being the most important due to the small orientation angle 

of ±10 degrees. 



Figure 11. Comparison between experimental data and both models for [+-10] s 
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Close examination of Fig. 1 1 reveals the superiority of the nonlinear model in predict- 
ing the composite behavior in the strain range of 0.05 to 0.10. However, the linear model 
seems to be more accurate in the strain range below 0.05; a feature left unexplainable due 
to the lack of sufficient experimental data. Nonetheless, the ability of the nonlinear model 
to capture the composite behavior was well demonstrated. 


5.0 Conclusion 

The micromechanical analysis of nylon reinforced elastomeric composites using 
MAC/GMC was greatly enhanced by accounting for the constituent’s material nonlinear 
behavior. The geometric nonlinearity due to small cord rotations during loading was 
assumed here to have a second order affect and consequently dropped from any further 
consideration. The user subroutine USRMAT within MAC/GMC, with analysis capabili- 
ties based on the generalized method of cells (GMC), was utilized to introduce the nonlin- 
ear constituent material behavior. Stress-strain behavior at the macro level was 
experimentally generated for single ply unidirectional and multi-ply symmetric laminate, 
continuously reinforced composites, composed of Nylon-66 embedded in a carbon black 
loaded rubbery matrix. Comparisons between the predicted macro composite behavior and 
experimental results, at both lamina and laminate levels were made with the agreement 
being very good provided the constituent’s material nonlinearity was included in the anal- 
ysis. Also, very good overall agreement with transverse strain measurements via Poisson’s 
ratio were obtained for all ply orientations when utilizing a nonlinear extension of the 
classical elasticity constitutive model. It is believed that the improvement realized at the 
macro level will also be observed at the micro level, a concept that will be further investi- 
gated in future works. 
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Appendix A 

Input file for the regular version of MAC/GMC using the unidirectional option 

Model 0 degree nylon 0.<e<0. 10 

* PRINT 

NPL=0 % 

*LOAD 

LCON=2 LOP=l LSS=1 % 

*MECH 

NPTW=2 TI=0., 1 . LO=0.,0. 10 % 

* MODEL 
MOD=l % 

* SOLVER 

NTF=1 NPTS=2 TIM=0.,1. STP=0.01 % 

*FIBER 

NFIBS=1 

NF=1 MF=6 MAT=6 IDB=n & 

EL=0.28236E6, 0.282363E5, 0.65 ,0.65, 0. 1E5,0.,0. % 

* MATRIX 
NMATX=1 

NM=1 MM=6 MAT=A IDB=n & 

EL=1.87E3, 1 ,87E3,0.49,0.49,0.690E3,0.,0. % 

*MRVE 

IDP=1 1 VF=0.247 RAD=0.0155 R=0.977374 % 

*CURVE 
NP=1 % 

* MACRO 
NT=1 

NC= 1 X=1 Y=7 NAM=nylonregularoutMAC % 

*END 
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A modified input file that utilizes the materials nonlinearity supplied by the user sub- 
routine USRMAT: 

Model 0 deg. hyten 0.<e<0. 10 
*PRINT 

NPL=0 % 

*LOAD 

LCON=2 LOP=l LSS=1 % 

*MECH 

NPTW=2 TI=0.,1. L0=0.,0.10 % 

*MODEL 

MOD=3 MATSYS=1 NLY=1 THK=0.03583 & 

CON=2 SYS=1 ANG=0. % 

*SOLVER 

NTF=1 NPTS=2 TIM=0.,1. STP=0.01 % 

*FIBER 

NFIBS=1 

NF=1 MS=1 MF=99 NPE=5 EL=0.28236E6,0.28236E5,0.65,0.65,0. 1F5& 

NPV=0 % 

^MATRIX 
NMATX= I 

NM=1 MS=1 MM=99 NPE=5 EL=1.87E3,1.87E3,0.49,0.49,0.690E3 & 

NPV=0 % 

*MRVE 
IDP= 1 1 

L=1 VF=0.2476 RAD=0.0075 R=0.977374 % 

*CURVE 
NP=1 % 

*MACRO 

NT=1 

NC=1 X=1 Y=7 NAM=modplyout % 

END 
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Appendix B 


Input file for the regular version of MAC/GMC using the laminate option 

Nylon +10/-10/+10/-10 laminate simulation 

* PRINT 

NPL=- 1 % 

*LOAD 

LCON=2 LOP=l LSS=1 % 

*MECH 

NPTW=2TI=0.,1. L0=0.,0.10 % 

*MODEL 

MOD=3 MATSYS=1 NLY=4 THK=0.055,0.055,0.055, 0.055 & 
CON=2,2,2,2 SYS=1, 1,1,1 ANG=+10.,-10.,+ 10.,-10. % 

*SOLVER 

NTF=1 NPTS=2 TIM=0.,1. STP=0.01 % 

*FIBER 

NFIBS=1 

NF=1 MS=1 MF=6 MAT=6 IDB=n & 

EL=0.28236E6, 0.28236E5 ,0.65, 0.65,0. 1E5,0.,0. % 

NPV=0 % 

* MATRIX 
NMATX= 1 

NM=1 MS=1 MM=6 MAT=D IDB=n & 

EL= 1 .87E3, 1 .87E3,0.49,0.49,0.690E3,0.,0. % 

*MRVE 

IDP=1 1,11,11,11 

L=1 VF=0.247 RAD=0.0155 R=0.977374 
L=2 VF=0.247 RAD=0.0155 R=0.977374 
L=3 VF=0.247 RAD=0.0155 R=0.977374 
L=4 VF=0.247 RAD=0.0155 R=0.977374 % 

*CURVE 
NP=1 % 

* MACRO 
NT=1 

NC= 1 X=1 Y=7 NAM=nylonlaminatecontrol % 

*END 
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The modified input file that utilizes the materials nonlinearity supplied by the user sub- 
routine USRMAT: 


Nylon +10/-10/+10/-10 laminate simulation 
*PRINT 
NPL=- 1 % 

*LOAD 

LCON=2 LOP=l LSS=1 % 

*MECH 

NPTW=2TI=0.,1. L0=0.,0.10 % 

*MODEL 

MOD=3 MATSYS=1 NLY=4 THK=0.055,0.055,0.055,0.055 & 

CON=2, 2,2,2 SYS=1, 1,1,1 ANG=+10.,-10.,+ 10.,-10. % 

*SOLVER 

NTF=1 NPTS=2 TIM=0.,1 . STP=0.01 % 

*FIBER 

NFIBS=1 

NF=1 MS=1 MF=99 NPE=5 EL=0.28236E6,0.28236E5,0.65 ,0.65,0. 1E5 & 
NPV=0 % 

*MATRIX 
NMATX= 1 

NM=1 MS=1 MM=99 NPE=5 EL=1 .87E3,1.87E3,0.49,0.49,0.690E3& 
NPV=0 % 

*MRVE 

IDP= 11,11,11,11 

L=1 VF=0.247 RAD=0.0155 R=0.977374 
L=2 VF=0.247 RAD=0.0155 R=0.977374 
L=3 VF=0.247 RAD=0.0155 R=0.977374 
L=4 VF=0.247 RAD=0.0155 R=0.977374 % 

*CURVE 
NP=1 % 

*MACRO 

NT=1 

NC=1 X= 1 Y=7 NAM=nylonlami 
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